Fracturing the optimal paths 
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Optimal paths play a fundamental role in numerous physical applications ranging from random 
polymers to brittle fracture, from the flow through porous media to information propagation. Here 
for the first time we explore the path that is activated once this optimal path fails and what 
happens when this new path also fails and so on, until the system is completely disconnected. In 
fact numerous applications can be found for this novel fracture problem. In the limit of strong 
disorder, our results show that all the cracks are located on a single self-similar connected line of 
fractal dimension Db ~ 1.22. For weak disorder, the number of cracks spreads all over the entire 
network before global connectivity is lost. Strikingly, the disconnecting path (backbone) is, however, 
completely independent on the disorder. 

PACS numbers: 62.20.mm, 64.60.ah, 05.45.Df 



CZ3 



I 

O 

o 



> ■ 
cn : 
vn 

p ; 

o ; 

On , 

o • 



X 



The identification and characterization of the optimal 
path in a disordered landscape represents an important 
problem in theoretical and computational physics as it 
can be intimately associated with many relevant scien- 
tific and technological applications 0, 0, H, 0, H, [1, 0, Hi 
including brittle fracture, random polymers and trans- 
port in porous media. It has been shown that optimal 
paths extracted from energy landscapes generated with 
weak disorder are self-affine and belong to the same uni- 
versality class of directed polymers Q . In the strong dis- 
order or ultrametric limit, on the other hand, numerical 
works U , l2| demonstrated the self-similar nature of 
the optimal path on two and three-dimensional lattices 
with fractal dimensions given by Df « 1.22 and 1.43, 
respectively. 

Optimal paths can be chosen by nature as those of min- 
imum energy as occurs for instance in electrical or fluid 
flow through random media or they can be chosen delib- 
erately in man made devices to reduce a cost function as 
for instance in internet or car traffic. Since these paths 
are heavily used they are more likely prone to failure 
either by overload, overheating or congestion. Two ques- 
tions that naturally arise are (?) how and when the system 
will eventually collapse and (zz) how the topology and in- 
homogeneity of this fracture affects its performance. It 
is the aim of the present Letter to provide a novel mod- 
elization for this complex problem that captures its essen- 
tial features and gives some insight about the statistical 
physics of these important questions. We flrst introduce 
a new model to generate a macroscopic fault line, namely 
the Optimal Path Crack (OPC) method, that is based on 
an iterative application of the Dijkstra algorithm [l3| to 
two-dimensional random landscapes. We then quantify 
the effect of disorder on the resulting crack topology to 
substantiate the relevance of the OPC method as a tool 
to understand the fracture of the random medium. 

Our substrate is a square lattice of size L with fixed 
boundary conditions at the top and bottom, and periodic 
boundary conditions in the transversal direction. Disor- 



der is introduced in our model by assigning to each site i 
an energy value ti ~ exp[/3(pi — 1)] [l3| , where pi is a ran- 
dom variable uniformly distributed in the interval [0, 1], 
and /? is a positive parameter that has the physical mean- 
ing of inverse temperature. It can be readily shown that 
this transformation is equivalent to choosing the values 
of e from a power-law distribution, namely P{^i) ^/^i, 
subjected to a maximum cutoff given by emax = ■ The 
existence of such a cutoff turns the hyperbolic distribu- 
tion normalizable for any finite value of (3. In this way, 
(3 represents the strength of disorder since it controls the 
broadness of the energy distribution. 

The energy of any path in the system is defined here 
as the sum of all the energies of its sites. In particular, 
the optimal path is the one among all paths connecting 
the bottom to the top of the lattice that has the smallest 
sum over all energies of its sites. This definition is simi- 
lar, but different from the definition of the shortest path 
in a network [l^. Without loss of generality, since we 
only consider positive pi values, the Dijkstra algorithm 
becomes a suitable tool for finding the optimal path [l^ . 

The OPC is formed as follows. Once the first optimal 
path connecting the bottom and the top of the lattice 
is determined, we search for its site having the highest 
energy which then becomes the first blocked site, i.e., that 
can no longer be part of any path. This is equivalent to 
impose an infinite energy to this site. Next, the optimal 
path is calculated among the remaining accessible sites 
of the lattice, from which the highest energy site is again 
removed, and so on. The blocked sites can be viewed 
as "micro-cracks" , and the process continues iteratively 
until the system is disrupted and we can no longer find 
any path connecting bottom to top. 

In Fig. la we show the resulting spatial distribution 
of blocked sites that constitutes the OPC in a typical 
random landscape generated under weak disorder condi- 
tions (/3 — 0.002). As depicted, the OPC structure has 
three basic elements. Besides the loopless backbone of 
the fracture (shown in blue) which effectively "breaks" 
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FIG. 1: (Color online) Typical realizations of the OPC model 
on a 512 x 512 lattice. In (a) the blocked sites generated un- 
der weak disorder conditions (/3 = 0.002) can be categorized 
in three distinct types, namely, the loopless backbone of the 
OPC (shown in blue) at which multiple dangling ends (shown 
in red) are attached, and isolated clusters (shown in black) 
spread all over the network. As shown in (b), under moder- 
ate disorder conditions (/3 = 6.0), while the OPC backbone 
is preserved, the number of dangling ends and isolated clus- 
ters is substantially reduced. Under conditions of very strong 
disorder (not shown), only the OPC backbone remains. 



the system in two, we can observe the presence of dan- 
gling ends (shown in red) as well as isolated clusters ho- 
mogeneously distributed over the entire network (shown 
in black). The situation becomes very different when we 
increase the value of the disorder parameter (3. As shown 
in Fig. lb, the amount of dangling ends and isolated clus- 
ters in a OPC generated under moderate disorder condi- 



tions (/3 = 6.0) becomes significantly smaller than in the 
case of weak disorder. By increasing further the value of 
/3, finally only the backbone remains. Interestingly this 
backbone is identical for all values of /3, while the en- 
tire set of blocked sites is highly dependent on the way 
disorder is introduced in the system. 
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FIG. 2: (Color online) Logarithmic dependence on size of 
the mass of all blocked sites Mt (circles) forming the OPC, 
the cluster of all sites in the fracture Mf (triangles), and the 
backbone mass M}, of the OPC that divides the system in 
two (stars), for /3 = 0.002. The system is considered to be in 
the weak disorder regime for this value of /3 and this range 
of system sizes. The three solid lines are the least-squares 
fits to the data of power-laws, Mt — L^* , Mf ~ L^i and 
Mb ~ with exponents Dt = 2.00±0.01, Df = 1.59±0.02, 
and Db — 1.22 ± 0.02, respectively. 

The behavior shown in Fig. 1 is somehow related to the 
problem of minimum path in disordered landscapes [l^ . 
In that case, the passage from weak to strong disorder in 
the energy distribution reveals a sharp crossover between 
self-affine and self-similar behaviors of the optimal path. 
In the strong disorder regime the energy of the minimum 
path is controlled by the energy of a single site. This sit- 
uation occurs when the distribution of energies can not 
be normalized, for instance, in the case of a power-law 
distribution, P{£i) ^ e~", for a < 1. The parameter 
(3 alone, however, does not determine the limit between 
weak and strong disorder, since this property also de- 
pends significantly on the system size. More precisely, 
if f3 is sufficiently high, or the lattice size is sufficiently 
small, the sampling of the distribution near the cutoff 
region is not so relevant. For any practical purpose, this 
network is considered to be in the strong disorder regime, 
resulting in a self-similar type of scaling for the minimum 
path. By increasing the network size, we may reach the 
point where one should expect to start sampling larger 
energy values that are beyond the cutoff of the distribu- 
tion. Above this scale, the system will return to the weak 
disorder regime, leading to a self-affine behavior for the 



3 



minimum path. In this way we should expect an abrupt 
transition from the weak disorder regime, at small values 
of 0, to the strong disorder regime, at large values of (3 
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In order to quantify macroscopically the effect of dis- 
order on the geometry of the OPC, we performed com- 
puter simulations for 1000 realizations of lattices with 
sizes varying in the range 32 < L < 512, and gener- 
ated with different values of the disorder parameter (3. 
In Fig. 2 we show for (3 = 0.002 that the average mass of 
the OPC backbone scales as Mf, ~ L^'' , with an exponent 
Db ~ 1.22 ± 0.02 [l^. Surprisingly, this exponent value 
is statistically identical to the fractal dimension previ- 
ously found for the optimal path line, but obtained under 
strong disorder l^, 11, l^j. It is also very close to that 
found "strands" in Invasion Percolation (1.22±0. 01 1C|), 
and paths on Minimum Spanning Trees (1.22 ±0.01 17|). 
In our case, however, it is important to note that the 
value of Df, reflects a highly non-local property of the 
system that is intrinsically associated with the iterative 
process involved in the OPC calculation. Also shown 
in Fig. 2 is the mass of the OPC fracture in weak disor- 
der, which consists of both the backbone and its dangling 
ends. While the crack itself grows as a power-law with 
size, Mf ~ L^f, with an exponent Df = 1.59 ±0.02, the 
total mass of blocked sites (crack and isolated clusters), 
however, is a constant fraction of the total mass of the 
system, i.e. Mt ~ L^. 
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FIG. 3: (Color online) Log-log plot of the masses Mt (circles), 
Mf (triangles) and Mb (stars) against system size L for /3 — 
6.0. For this intermediate value of /3, we can clearly identify 
the crossover from strong to weak disorder depending on the 
system scale. For small lattices, the system operates under 
strong disorder conditions. As a consequence, most of the 
blocked sites lie on the fracture path and the three masses Mt , 
Mf and Mb are identical. As the system size increases, one 
reaches the weak disorder regime, and the three curves split 
apart. At larger scales, we should recover the same power-law 
behaviours found for sufficiently low values of (3, as shown in 
Fig. 2. 



In Fig. 3, the results obtained for P = 6.0 clearly indi- 
cate the transition from strong to weak disorder regimes 
by systematically increasing the size of the system. As 
already mentioned, the stronger the disorder in the sys- 
tem (low L or high /?), the smaller is the number of final 
blocked sites that also become more and more localized 
in a singly-connected crack line. Precisely, in the limit of 
very strong disorder, we obtain that only the OPC back- 
bone mass Mf, remains, i.e. Alt — > Mi, and Mf M^, 
scaling in the same way as in the weak disorder limit, 
namely Mb ~ L^" , with Db = 1.22 ± 0.02. As shown in 
Fig. 1 , the backbone is indeed invariant under the change 
of (3. By increasing the value of L, we observe a gradual 
departure of and Alf from this behavior to their re- 
spective scaling laws in weak disorder (high L or low f3) , 
as displayed in Fig. 2. 
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FIG. 4: (Color online) Transition from weak to strong disor- 
der in the OPC model. In these curves we plot the density 
of all blocked sites pt at the end of the process when the top 
of the system disconnects from the bottom. These curves 
present three distinct regimes. For weak disorder (3 < 1 one 
sees a saturation in the density in a value around pt ~ 0.22. 
For intermediate values one sees the density decreasing with 
the disorder parameter /3, however there is no dependence on 
the system size. For larger values of (3 we see the density 
saturating again, but now to a value that depends on the sys- 
tem size L. This last regime is the one that characterizes the 
strong disorder regime. The saturated density in this regime 
depends on the system size as L^''~^, where Db ~ 1.22 is the 
fractal dimension of the OPC backbone. For each system size, 
we observe the transition to the strong disorder regime at a 
different value of the parameter f3x ■ From the collapse shown 
in the inset, a non-trivial dependence is revealed between the 
onset of the transition and the system size, /3x ~ L°'^^. 

The transition from weak to strong disorder is better 
illustrated by the results depicted in the inset of Fig. 4. 
There we plot the density of all blocked sites pt at the 
end of the process as a function of the parameter /?. The 
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curves exhibit three different regimes, depending on the 
value of p. For small values, (3 < 1, the density saturates 
to a fixed value. For larger values, the density decays as 
a power-law pt with an exponent 9 « 4/3. Both 

regimes, saturation and power-law decay, are still in weak 
disorder. For finite lattices, the curves present another 
crossover to a minimum density that is now dependent 
on the system size. This second crossover, /3x , which 
indicates the transition to strong disorder, should depend 
on the system size in such way that an infinite system 
is in weak disorder for any finite value of (3. For large 
enough values of /3, in the strong disorder regime, the 
density reaches a minimum value when all the blocked 
sites lie on the fracture dividing the system. As shown 
in Fig. 2, the mass of this fracture scales as L^'' . At the 
onset of the transition, this power-law behavior crosses 
over to a scale dependent value. Thus, ~ L^''~^, 
or /5x i", with a = (2 - Db)/0. The collapse of the 
results for intermediate and large values of P obtained 
using Df) = 1.22 and a — 0.59 shown in the main plot of 
Fig. 4 is consistent with this analysis. 

Concluding, we discover that for all disorders the line 
along which all minimum energy paths fracture is fractal 
of dimension 1.22 in 2d. The role of disorder and system 
size can be fully cast in a crossover scaling law for the 
total number of blocked sites. Our model poses new chal- 
lenges also from the theoretical point of view since the 
numerical resemblance of our fractal to the one of domain 
walls and optimal paths in the strong disorder limit [l(j | 
seems to hint towards some deeper relation. It could be 
certainly interesting to study also the influence of the 
dimension of the system and of the substrate topology 
on our model by a generalization to three dimensional 
systems or complex network. 
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